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K.  Y.  Kim 

E.R.L.,  University  of  California,  Berkeley, CA.  94720 


ABSTRACT 

JWe-  present^wa/yr/c  solutions  for  asymmetric  double  layers  which  satisfy 
the  time  stationary  Vlasov-Poisson  system  and  which  require  the  double¬ 
valuedness  of  Sagdeev  potential  as  a  function  of  physical  potential:  it  is  pointed 
out  that  any  distribution  function  having  an  analytic  density  representation  as  a 
polynomial  power  series  of  potential  can  never  satisfy  the  asymmetric  double 
layer  boundary  conditions.  Considering  K-dV  like  equation,  it  is  found  that 
there  is  some  relationship  between  the  speed  of  asymmetric  double  layer  and 
the  degree  of  asymmetry.  c~. 

A  monotonic  double  layer  is  simply  an  isolated  pair  of  oppositely  charged  sheets  which 
result  in  a  narrow  region  of  abrupt  potential  jump  of  some  amplitude  —  <fr;  well  outside  of 
this  localized  jump,  the  potential  is  effectively  uniform1'20.  An  asymmetric  double  layer  refers 
actually  to  a  localized  region  of  three  sheets  of  alternating  charge  sign,  having  an  asymmetric 
potential  jump  profile  and  thus  including  sub  regions  of  oppositely  directed  asymmetric  electric 
fields.  Asymmetric  double  layers  were  first  reported  in  recent  numerical  simulations  and  in 
space  observations  on  auroral  field  lines9"11;  subsequent  theoretical  studies  have  numerically 
examined  the  time  stationary  Vlasov-Poisson  system7-8-12.  It  has  been  suggested  that  small 
amplitude  double  layers  may  account  for  a  large  portion  of  the  total  potential  along  auroral  field 
lines  and  may  also  explain  the  fine  structure  of  auroral  kilometric  radiation:  recent  satellite 
measurements11  are  especially  consistent  with  the  asymmetric  double  layer,  its  potential  depres¬ 
sion  (or  the  potential  hump)  at  the  low  potential  side  (or  high  potential  side).  Recent  studies  of 
the  thermal  barrier  in  tandem  mirror  devices  have  also  found  that  there  exist  states  with  abrupt 
potential  depressions  as  a  result  of  forced  changes  in  the  distribution  functions7.  The  negative 
potential  depression  is  thought  to  play  a  crucial  role  in  the  formation  of  double  layers  and  to  be 
responsible  both  for  current  disruptions  (by  reflecting  the  electrons)  and  also  for  high  fre¬ 
quency  noise  excitation  behind  the  double  layer  (by  a  two  stream  instability  involving  electrons 
that  pass  the  negative  potential  peak8"11 ). 

Although  there  have  been  many  theoretical  and  experimental  investigations  of  double 
layers,  recent  theoretical  work  has  been  limited  to  only  numerical  evaluations  of  the  Vlasov- 
Poisson  system  (or  of  the  fluid  system)  mainly  because  of  the  highly  nonlinear  properties  of 
asymmetric  double  layers.  One  such  recent  numerical  calculation12  suggested  that  there  may  be 
a  low  amplitude  limit  for  the  monotonic  double  layer.  In  a  previous  paper20,  we  gave  analytic 
evidence  for  the  existence  of  small  amplitude  monotonic  double  layers,  which  are  the  analytic 
extensions  of  the  well  known  electron  solitary  hole  and  ion  acoustic  solitary  hole21-22.  In  this 
report,  we  present  the  self-consistent  analytic  solution  for  asymmetric  double  layers,  which 
satisfy  the  time  stationary  Vlasov-Poisson  system.  We  present  a  derivation  of  a  K-dV  tike  equa¬ 
tion  describing  a  moving  asymmetric  double  layer.  It  is  shown  that  asymmetric  double  layer  has 
a  some  relationship  between  the  asymmetry  parameter  and  the  speed  of  double  layer. 

To  describe  propagation  of  an  electrostatic  double  layer,  we  use  a  Vlasov-Poisson  system 
that  has  been  Galilean-transformed  to  the  wave  frame  (where  the  wave  is  time  stationary).  In 
this  frame,  we  can  express  the  time  stationary  solution  of  Vlasov  equation  as  any  function  of 
the  constants  of  motion:  (i)  particle  total  energy  and  (ii)  the  sign  of  the  velocity  of  the 
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ua trapped  particles23.  Besides  the  above  usual  constants  of  motion,  it  is  important  to  realized 
that  a  third  constant  of  motion  exists  for  the  reflected  particles:  sgn  (  x  -  xm  )  where  xm 
represents  the  position  of  potential  minimum  (or  maximum)  for  the  negatively  charged  parti¬ 
cles  (or  the.  positively  charged  particles).  It  turns  out  that  this  final  constant  of  motion  plays  an 
important  role  in  order  to  construct  the  asymmetric  double  layers.  With  the  above  distribution 
function(any  function  of  the  three  constantsof  motion)  ,  we  can  develop  a  general  theory  of 
double  layers  by  expanding  the  corresponding  density  in  half-integer  powers  of  ±4-  However, 
in  this  paper  we  would  like  to  give  a  physically  more  transparent  description  of  our  double  layer 
problem,  so  we  will  proceed  by  using  specific  distribution  functions. 


In  order  to  describe  the  asymmetric  double  layer  accompanying  a  potential  depression  at 
the  low  potential  side,  we  assume  a  Maxwell- Boltzmann  ion  distribution, 
h.'Jt  —  »2  +  2*  > 

2  ,  and  we  consider  the  following  electron  distribution  function  . 


f,  -  (2ir)-*  (  expMOgn  v  «*-  vd  )2  }  ©(a) 

+  ( 1  +  /  &»(  x  -  xm  ) )  exp{-ft(v/  +0  «)  }  ©(-«)  ] 


(1) 


where  r  -  TJT,  and  <  -  v2  -  24  for  0  <  4  <  4  .  Here  the  electron  velocity,  the  wave 
potential  and  the  spatial  coordinates  are  normalized  to  the  electron  thermal  velocity  (7’,/m,)  *  , 
the  electron  temperature  Tje  and  the  electron  Debye  length  A,  -  ( Tt/^rrn0e2)  *,  respectively; 
vd  represents  the  electrons  drift  velocity.  The  ion  velocity  has  been  normalized  by  the  ion 
acoustic  velocity.  Here  H,  and  /  are  constants  and  0  represents  an  inverse  temperature.  Here 
-I  <  /  <  1  in  order  to  obtain  positive  trapped  particle  distribution  function. 

With  these  distribution  functions  /,  and  /,,  the  Poisson  equation  may  be  written  in  the 
small  amplitude  limit  by  introducing  a  Sagdeev  potential  K(4)  : 
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-1-  {  higher  order  terms  } 


where  Ed  is  the  electron  drift  energy  and  Z*,(x)  represents  the  real  part  of  the  derivative  of  the 
complex  plasma  dispersion  function(see  Fig.  1.  (a)).  The  double  layer  solution  may  be  found 
by  considering  the  following  nonlinear  eigenvalue  conditions: 

Charge  neutrality  at  x  -  ±«>  requires  that  the  rhs  of  Eq.(3)  should  vanish  at  the  boun¬ 
daries  4  -  <fri ,  4* 

*  Existence  of  the  asymmetric  double  layer  requires  that  the  Sagdeev  potential  be  identically 
zero  at  4  m  0  ,  4| ,  4:>  *o  that  the  electric  field  equals  zero  at  those  values  of  potential. 

•  An  additional  condition  on  the  Sagdeev  potentiallsee  Fig.  2.  (a)}  requires  that  K(4)  <  0 
for  0  <  *  <  4i ,  0  <  *  <  *2. 

Here  it  is  important  to  note  that  the  Sagdeev  potential  should  be  a  double  valued  function 
of  *  for  0  <  *  <  min  (4 1 ,  #2)  in  order  to  yield  the  asymmetric  double  layer  solution,  while 
the  existence  of  the  monotonic  double  layer  and  the  solitary  hole  does  not  require  the  double 
valuedness  of  Sagdeev  potential.  Furthermore,  it  should  be  noted  that  from  the  double  valued- 
ness  of  Sagdeev  potentiaKas  a  function  of  4)  for  the  existence  of  the  asymmetric  double  layer. 


we  can  make  a  following  statement:  any  distribution  function  having  analytic  density  represen¬ 
tation  as  any  polynomial  power  series  of  potential  («*J*  /  dv  -  £  C„  (  ±<t>  )  with 

— «•  It 

C  *■  constant )  can  never  satisfy  the  asymmetric  double  layer  boundary  conditions.  In  our 
^  double-valuedness  of  Sagdeev  potential  is  guaranted  by  the  use  of  the  constant  of  motion 
for  reflected  particles. 
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FtO.  1.  U)  The  real  put  of  derivative  of  the  complex  plasma  dispersion  function:  —Z'r(x).  (b)  The  Sagdeev  poten¬ 
tial  for  the  asymmetric  double  layer. 

By  solving  Eqs.(2)  and  (3)  subject  to  these  conditions,  we  get  the  following  double  layer 
solution: 

4  «■  £{  o,  +  tanh  ±|k,  |  x  P 

where  for  convenience  consider  only  >  <l»i  and  we  have  defined 

(•*/*!  +  yfih)2  Vfo  ~  ■>  v * 

° 5  •v/'fri  +  ’ 


*- 


-  -H5L  >  0 


(4) 


(5) 


4  ’  V*l  +  V*2  ’  "  6 

Here  *,  a„  k,  are  related  to  our  system  parameters  as  follows: 

l-«,  -  v  * 

f  e  E*  *  ±  “  V5:)  *  * 

{  «,r  -  iz’,(v^)  )  “  y(*i  +  *2  “  *  . 

»-( ju-ft'1  ,£‘ > +  Ifr  > 

For  example  if  we  set  -  0  (a,  -  0)  in  Eq.(4),  then  we  would  obtain  our  previous 
monotonic  electron  double  layer  solution20,  which  is  the  nonlinear  version  of  slow  electron 
acoustic  wave  in  the  limit  r— 0,  #2  —0.  Similarly,  the  condition  #•  -  *:  (a,  -  0)  corresponds 
to  the  solitary  ion  hole  solution.  More  generally  the  condition  M<1,  the  above  solution 
Eq.(4)  describes  an  asymmetric  double  layer,  with  a  potential  depression  at  the  low  potential 
side.  Therefore,  we  may  consider  a ,  to  be  the  asymmetry  parameter.  From  Eq.(6),  we  see  that 


(6) 
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(9) 

(10) 
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the  existence  of  the  asymmetric  double  layer,  with  a  potential  depression  at  the  low  potential 
side,  requires  a,  <  l  in  order  to  have  a  positive  curvature  at  <£  -  0. 

Having  considered  the  asymmetric  double  layer  with  a  potential  depression  at  the  low 
potential  side,  we  now  turn  to  the  problem  of  an  asymmetric  double  layer  having  a  potential 
hump  at  the  high  potential  side.  In  order  to  describe  this  class  of  double  layer,  we  assume  a 

-  -f-Crv2  -  2*) 

Maxwell-Boltzmann  electron  distribution  /,  -  ntun)~*e  2r  and  we  consider  the  fol¬ 

lowing  ion  distribution24: 

/,  -  (2»)-*|  exp{—  sgn  »  t*  +  v0  )2}  9 («,) 

+  (  l  +  A  Sg/i  (  x  -  xm  )  )  exp  {— #ae,}  8(— «,)  j  (11) 

where  «,  -  vJ  +  2r<p  with  H  <  <t>  <  0  ,  r  -  TJ  T,  and  ion  velocity  has  been  normalized 

to  the  ion  thermal  velocity  ( TJm ,)*  .  Here  v0  and  £0  represent  the  ion  drift  velocity  and  the 
ion  drift  energy,  respectively.  Here  H,  and  h  are  constants  and  a  represents  an  inverse  tem¬ 
perature.  Here  —  1  <  h  <  1  in  order  to  have  positive  trapped  particle  distribution  function. 
With  the  above  distribution  functions,  the  small  amplitude  limit  Poisson  equation  becomes 

djop  "  («,-  l  )  +  -  h  -  xm  H~d)1/J  +  {  h,  -  J  Z'rlJTo  )  )<* 

(  1+A  sgn(x  -  x^)  )  a  -l)(H)3/: 

+  t2{  -  e~£#(j 6  +  -|^y)  l*2  +  (  higher  order  terms  )  (12) 

Charge  neutrality  at  x  -  ±°°  requires  that  K*(^)  -  0  at  4  »  — ,  Hj.  The  conditions  for 

zero  electric  field  lead  to  K(^)  -  0  at  4  —  0  ,  Hi  ,  Hj.  Again  solving  Eq.(12)  and  the 
above  nonlinear  dispersion  relations  together  with  the  requirement  K(6)  <0  {see  Fig.  2.  (b)} 
for  <  4  <  0  ,  Hi  <  4  <  0,  we  obtain  the  following  asymmetric  double  layer  solution 
with  a  potential  hump  at  the  high  potential  side(|a, |  <  1): 

4  «■  H{  a,  +  tanh  ±  |k,|  x  }2  (13) 

where  for  convenience  we  have  defined 

r  (y/+\  +  V*j)2  _2_  /u\ 

* - 4  'a'  7XT7Z'*'-  6-  (14) 

Here  #t,  and  k,  are  related  to  our  system  parameters  as  follows: 

-  1  -  ,  05) 

-  ±  .  06) 

{  A,  -  |tZ*,  (V^)  )  -  -y <*,  +*:-  ,  (17) 

~ *  *#Aa  -  T  y  (V*7  -  ,  08) 

M.|A_^-«.(X  +  20,I.  (i» 

For  example  if  we  set  #,  •  0  (a,  —  1)  in  Eq.(10),  then  we  would  recover  our  previous 
ion  acoustic  aoootooic  double  layer  solution20,  which  is  the  nonlinear  version  of  the  slow  ion 
acoustic  wave.  Similarly  the  condition  #]  -  (a,  -  0)  would  return  the  electron  solitary  hole 
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solution.  From  Eq .(10),  we  note  that  the  existence  of  the  asymmetric  double  layer  with  a 
potential  hump  at  the  high  potential  side,  requires  H,  <  1  in  order  to  have  a  negative  curvature 

at  4  •  0. 

Having  obtained  the  analytic  solutions  for  the  time  stationary  double  layers,  we  can  now 
present  a  derivation  of  the  K-dV  like  equation,  which  describes  the  one  dimensional  asymptotic 
behavior  of  the  asymmetric  double  layers,  having  a  potential  depression  at  the  low  potential 
side,  of  small  but  finite  amplitude. 

To  describe  a  collisionless  plasma  of  cold  ions  and  warm  electrons,  we  consider  the  fol¬ 
lowing  set  of  equations: 

n,  4-  (nv)x  -  0  ,  v,  +  w,  +  -0  ,  -  n,  -  n  ,  (20) 


where  the  density,  velocity,  potential  and  spatial  coordinate  are  normalized  to  the  unperturbed 
density  n<y,  ion  acoustic  velocity  ( Te/m,)ui ,  the  electron  temperature  Tje  and  the  electron 
Debye  length,  respectively.  Assuming  the  electrons  to  be  in  a  quasi-equilibrium  with  the  low- 
frequency  ion  acoustic  wave,  we  may  expand  the  electron  density  as  before 
it,  —  1  +  c0  sgn  (  x  —  xm  )  ($)I/J  +  C\4  +  cj  sgn  (  x  -  xm  )  (6)}/:  +  c302  + .  By  intro¬ 

ducing  the  Gardner-Morikawa  coordinate  transformation  S  —  $1/2  (x—t)  and  r  -  bylt  ,  then 
expanding  n,v,d  in  powers  of  small  parameter2*  8,  and  then  assuming  c0  and  c2  to  be  order  of 
8V2  and  8l/2,  respectively,  we  obtain  the  following  K-dV  like  equation  by  using  the  boundary 
conditions  for  an  asymmetric  double  layer  with  a  potential  depression  at  the  low  potential  side 
as  before: 


0- 


1  A  4.  — 1  ■ 
*«2*'  2*«2 


(21) 


-  ±  6 a,  (1-a,  *)  sgn  (x  -  x„ )  (-£■) 


1/2 


+  {4(3a,2— 1)  -  -£)£  *  10*  *"<*  "  ^  +  3(^): 

K  ih  ilr  dr 


where  M  represents  the  velocity  of  the  double  layer,  a2  - 


* 

(Af-1) 


and  >  0.  It  should 


4(3a,2  —  1) 

be  noted  that  we  have  used  our  boundary  condition  for  a  moving  asymmetric  double  layer  so 
that  we  can  extract  some  useful  physics.  Here  we  have  the  same  definition  of  i  and  a,  as  in 
Eq.(5).  The  evolution  equation  for  the  asymmetric  moving  double  layer  with  a  potential  hump 
at  the  high  potential  side  can  also  be  obtained  similarly  by  usinjj  the  corresponding  asymmetric 
double  layer  boundary  conditions;  it  is  given  by  simply  letting  <1/  —  — i  <  0  in  Eq.(21),  using 
the  definition  of  £  as  in  Eq.(14).  The  corresponding  moving  asymmetric  double  layer  solution, 
with  a  potential  depression  at  the  low  potential  side,  of  Eq.(21)  is  given  by 

2 


6(x ,/)  -  £{  a,  +  tanh  ± 


V 


Af-1 


4(3a,2  -  1) 


( x-Aft )  j 


(22) 


where  xm  is  given  by  the  equation  ^  (  xm  —  x  -  Mt ,  0  )  —  0.  Here  it  should  be  noted  that 
M  >  1  for  3a,2  -1  >  0  and  that  M  <  1  for  3a,2  -1  <  0  . 

In  conclusion,  we  have  obtained  two  different  asymmetric  double  layer  analytic  solutions: 
one  has  a  potential  hump  at  the  high  potential  side,  the  other  has  a  potential  depression  at  the 
low  potential  side.  By  considering  the  double-valuedness  properties  of  the  Sagdeev  potential 
(required  for  the  existence  of  an  asymmetric  double  layer),  we  have  proven  the  following  fact: 
any  distribution  function  having  analytic  density  representation  as  any  polynomial  power  series 

*9 

of  potential  (for  example  n  -  J  /  dv  -  £  C„  (  ±<8  )"  with  C„  -  constant )  can  never 

— m  * 

satisfy  the  asymmetric  double  layer  boundary  conditions.  This  shows  the  importance  of  using 
the  third  constant  of  motion  for  reflected  particles,  in  order  to  satisfy  the  double-valuedness  of 
Sagdeev  potential  for  the  asymmetric  double  layers.  We  have  also  given  a  derivation  of  the  K- 
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dV  like  equation,  which  describes  the  asymmetric  moving  double  layer  with  a  potential  depres¬ 
sion  at  the  low  potential  side26.  Thus  we  have  found  that  there  is  some  relationship  between  the 
speed  of  asymmetric  double  layer  and  the  degree  of  asymmetry  (a,). 

I  would  like  to  thank  to  Prof.  C.  K.  Birdsall,  Dr.  T.  Crystal,  Dr.  W.  Lotko,  Mr.  V.  Tho¬ 
mas  and  Mr.  N.  Otani  for  useful  discussions.  This  work  was  performed  under  the  auspices  of 
DOE  at  ERL,  University  of  California,  Berkeley. 
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